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This study is concerned with the diffusion of a passive scalar 8(r, t) advected by general n-dimensional shear 
flows u — u{y, z, ■ ■ ■ , t)x having finite mean-square velocity gradients. The unidirectionality of the incom- 
pressible flows conserves the stream-wise scalar gradient, ^xQ^ allowing only the cross-stream components 
to be amplified by shearing effects. This amplification is relatively weak because an important contributing 
factor, 9x6, is conserved, effectively rendering a slow diffusion process. It is found that the decay of the 
scalar variance (9^) satisfies d{&'^)/dt > —Ck^^^, where C > is a constant, depending on the fluid ve- 
locity gradients and initial distribution of 0, and k is the molecular diffusivity. This result generalizes to 
axisymmetric flows on the plane and on the sphere having finite mean-square angular velocity gradients. 

The transport of a diffusive scalar Q{r,t) by incompressible fluid flows u{r,t) is governed by 

{dt+u-V - KA)e = 0, (I) 
V-M = 0, 

where k is the molecular diffusivity. This system has been an active research subject^" for its application to 
problems in geophysical, environmental, and industrial context. Various types of flows have been considered, 
ranging from single-scale^'^ to fully turbulent flows. ^'^ A problem of genuine interest has been the decay 
behaviour of a finite initial scalar distribution in the limit of small k, where the question is whether the 
decay rate d{Q^)/dt of the scalar variance (9^) remains nonzero as k — > 0. In two dimensions, recent 
results^'^ have suggested a negative answer to this question. However, one of these results^ was derived for 
finite time only, while the other result^ relied on an assumption of power-law spectra of (9^). 

For shear flows, the advection-diffusion problem is more tractable, even in higher dimensions, making it 
possible for the above answer to be derived rigorously, as is presently shown. Moreover, in this case, the decay 
rate approaches zero relatively rapidly as k — > 0. The main feature in shear flows, not shared by turbulent 
flows, that facilitates this derivation is that shear flows conserve the stream-wise scalar gradient, say dxQ, 
allowing only the cross-stream components to be amplified. The amplification of these components is non- 
accelerated since an important contributing factor, dxQ, is conserved. More quantitatively, shearing effects 
can give rise to a quadratic growth of the mean-square scalar gradients (|V9p), as opposed to exponential 
growth, presumably occurring for turbulent flows. This slow growth of (|V9p) is overwhelmingly suppressed 
by diffusive effects, even in the limit k — > 0. As a consequence, given a smooth initial scalar distribution, 
the maximally achievable value of (|V9p) grows relatively slowly as k is decreased. This results in a rapid 
approach of the decay rate K(|V9p) to zero as k 0. 

This Brief Communication derives an upper bound for the decay rate K(|V9p} for general n-dimensional 
shear flows u = u(y, z, • • • , t)x. For simplicity, periodic boundary conditions are considered, and the flows 
are assumed to have finite velocity gradients or just finite mean-square velocity gradients. The initial scalar 
distribution is assumed to have finite mean-square gradients. The derived upper bound for K(|V9p) is found 
to scale as k.^^^. This analytic result is valid uniformly in time and generalizes to axisymmetric flows on the 
plane and on the sphere having integrable differential angular velocity. 

For u = u{y, z, - ■ ■ , t)x, equation (1) becomes 

{dt+udx-KA)e^O. (2) 

As u{y, z, • • • , t) is independent of the flow direction x, the classes of Fourier modes 9(fc, t) having a common 
kx are dynamically decoupled.^ This is directly connected to the fact that the derivatives 9™9, for m = 
1, 2, 3, ■ • ■, are materially conserved, as can be seen from the governing equation for 9™ 9, 

{dt + udx-KA)d'x''Q^Q. (3) 

This advection-diffusion equation is the same as that for 9; hence, all conservation laws for 9 also apply 
to 9™9. In particular, the suprcmum ||9a:9||^ and the mean-square stream-wise scalar gradient (|92;9p). 
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both being used in the subsequent calculations, are conserved by the advcction term. Under diffusive effects, 
these quantities decay in time, and hence are bounded by their initial values. 

The approach of Tran^^ (see also Tran and DritscheF) for estimating the enstrophy dissipation in two- 
dimensional turbulence is now applied to the present case. As the evolution of dxQ is trivial, one can drop 
this component from consideration. However, there appears to be no gain for so doing. Hence, for clarity, 
all the subsequent calculations make no separation of dxQ from V6. The evolution equation for V6 is 

{dt + udx- kA) Ve = -dxQ Vm. (4) 

The 'forcing' term on the right-hand side of (4) depends on the velocity gradient Vu and on the decaying 
stream- wise gradient 9x0, but does not depend on the cross-stream gradients. This means that each cross- 
stream component is amplified independently and that (|V0p) cannot grow more rapidly than quadratic in 
time. It is also worth mentioning that since d^Q can only decay, an initial scalar distribution homogeneous 
in the flow direction, i.e. d^Q = 0, decays purely diffusively because the right-hand side of (4) is identically 
zero for all t > Q. In other words, shearing effects alone cannot amplify any cross-stream scalar gradients 
in the absence of 9x0- This is obvious from physical point of view. The evolution equation for (|V6p) is 
obtained by multiplying (4) by V© and taking the spatial average of the resulting equation, 

\jtW) = -{d.Q . ve) - /^(lAep), (5) 

where the advection term identically vanishes. 

Two estimates of the triple-product term in (5) are obtained by using the Cauchy-Schwarz inequality: 

Either equation of (6) can be used to deduce an upper bound for k(1V0|^). depending on which estimate 
of the triple-product term is more optimal. For the first equation, applying the Cauchy-Schwarz inequality 

(|Aep) > (|ve|2)V(e2) yields 

^|(|veH < {\Vu\^Y'^ \\dM^ (ive^v^ - ^^^^^ (7) 

The first term in the brackets on the right-hand side of (7) depends on (jVup) (which is equivalent to 
(|V X itp) in this case) and on the decaying quantities ||9a;0||(^ and (0^). Given a bounded (|Vup) and 
bounded initial ||5a:6||oo and (O^), this term remains bounded. It follows that for f > 0, 

K(|ven3/2 < c, (8) 

provided that it holds for t = Q. Here, c is an upper bound for (|Vwp)^/^ ll^xQUoo (Q^)- From (8) one can 
readily deduce that 

where C is an upper bound for 2(|Vup)i/3 \\dx@tL^ (62)2/3. The decay of (G^) then satisfies 

^(62) = -2K(|Ve|2) > -Ck}/^. (10) 

The same bound but with a different constant C can be obtained by manipulating the second equation of 
(6) along these lines. In this case, C denotes an upper boimd for 2 ||Vu||^^ (|5xe|2)i/3(e2p/3^ instead 
of the requirement (|Vu|2) < cxd, the slightly more stringent condition ||Vu||j^ < oo is assumed. 

Equation (10) implies a slow decay of (92) in the limit of small /t. In other words, shear flows are 
rather ineffective mixers. The slow decay of (62) has a bearing on its exponential decay rate, which is of 
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particular interest and has been widely studied. For handling exponential decay behaviour, equation (10) 
can be rewritten in the more convenient form, 

^ '^:{e')>-^. (11) 



(92) (92) 

Suppose that at any instance in time, the decay of (92) is approximated by a pure exponential decay at the 
rate of A, then by (11), A satisfies 

A < m 

This equation gives an explicit upper bound for A in terms of k. 

During the period (or periods) of scalar gradient growth, i.e. d{\VQ\^)/dt > 0, the constant C in (12) 
need not be an upper bound for the quantity in question, but can be its instantaneous value (cf. (6) and the 
subsequent calculations). Therefore, in such a period, A satisfies 

A - < / 2(|V^.|2)V3 ||a^e||^3 ^02^-1/3,1/3^ 

(92) - \2||V^^f^/^(|9,9|2)l/3(e2)-l/3,l/3, ^''') 

For the period (or periods) of scalar gradient decay, i.e. d{\VQ\'^)/dt < 0, the validity of (13) could become 
questionable only if (|V9|2) has decayed more slowly than (92), i.e. the small scales have decayed relatively 
more slowly than the large scales. This condition requires that the production of cross-stream scalar gradients 
remain considerable throughout their decay. In any case, from (12) one can conclude that in the limit of 
small K, A could become sizable only when (92) oc k^/"^, i.e. when the scalar distribution has become 
virtually homogeneous. Furthermore, the remaining fraction of (92) would then have been 'cascading' to 
a dissipation wavcnumber, say kd, satisfying ka oc k^^^^. It is notable that for a simple shear flow in two 
dimensions satisfying all the required conditions, Vanneste and Byatt-Smith^2 argue that fast decay of the 
scalar is possible, in the sense that A remains nonzero in the limit k ^ 0. The present result implies that 
the asymptotic regime for considering this possibility is (92) oc k^^^. 

The; presciiit result readily generalizes to axisymmetric flows on the plane (vortical flows) and on the 
sphere (zonal flows), provided that these flows do not become singular, in the sense to be described in due 
course. For these cases, the respective advection-diffusion equations in the plane polar and spherical polar 
coordinates are 

{dt + 0(r, t) 84, - K{r-^ dr{rdr) + r'^ a|)) 9 = (14) 

and 

[dt + 0(6»,i) 84, - K{sm-'^ e 80{sm.e de) + sin"^ 6»a|)) 9 = 0. (15) 

Here ri(-,t) denotes the fluid angular velocity, and all other notations are standard. The sphere radius has 
been set to unity. For (14) zero boundary conditions are imposed. Similar to the previous case, the stream- 
wise scalar gradient in each case (9,^9/r for vortical flows and 9,^9/ sin ^ for zonal flows) is conserved by 
the flows. However, these arc not known to decay under diffusive effects. Instead, the decaying quantities in 
these cases are the derivatives 9^*9, which are governed by the same equation as 9 in each case. Now it is 
straightforward to perform the above calculations in curvilinear coordinates. For vortical flows, one obtains 

ld_ 
2di 

< {\dre\y/'(\drn\y/' 119^911^ - ^(|A9|2) (16) 

(r^i/2 

- (92) 
Similarly for zonal flows, one obtains 
1 d 



oiidvei") = (Aend^Q) - k{\a&\'') = -(a,9a,fia^9) - K{\Aef) 



({\drn\y^' M|^(92)-K(|V9|V/2). 



^^(|V9|2) = (A9na^9) - «(|A9|2) = -{dee8end^e)-K{\Ae\^) 

< {\dee\y/'{\don\y/' ||9^9|U - ^^(|A9n (17) 

< ^^^^sr (^"^^""'^'^' 1'^^®"- - «(iv9p)3/^) . 
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Prom (16) and (17), one can recover the result d{&'^)/dt > —Ck^^^ derived earlier for the periodic case. The 
constant C is an upper bound for 2(|9^f^|2)i/3 ||a^e||^-'^ (92)2/3 and 2{\dgn\Y^^ P^©!!^^ (92)2/3 for the 
vortical and zonal cases, respectively. Here ||909||j^ plays the role of ||92;9||j^ and (l^r-i^P) and (jSerip) play 
the role of (|Vm|2). The condition for the flows is the integrability of the differential rotation (l^r-i^P) < oo 
and (l^gOp) < oc. rather than the integrability of the velocity gradients. Note that for the case of vortical 
flows, Bajer, Bassom, and Gilbert^'^ have found by a different method that the decay rate k(|V9|2) scales 
as fc-^/^ (also see Rhines and Young^'*). 

In conclusion, this Brief Communication has derived a rigorous upper bound for the decay rate of the 
variance (92) of a passive scalar 9 in general n-dimensional shear flows. The flows are assumed to have flnite 
velocity gradients or just finite mean-square velocity gradients, and the initial scalar distribution is assumed 
to be smooth. This upper bound is valid uniformly in time and scales as k^/'^, where k is the diffusivity. 
This implies that in the limit of small diffusivity, the diffusion of a passive scalar in shear ffows is slow: shear 
flows are rather poor mixers. The reason is that shear flows conserve the stream- wise scalar gradient and 
amplify the cross-stream scalar gradients relatively weakly. This residt generalizes to axisymmetric flows on 
the plane and on the sphere having finite mean-square angular velocity gradients. 
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Bassom, and Gilbert^^ and the extension from the Cartesian to curvilinear cases. He would also like to 
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Pierrehumbert and an anonymous referee for comments, which were helpful in improving this manuscript. 
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